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Abstract 

We examine the spectrum of small perturbations around global and local (gauge) 
abelian vortices, using simple numerical matrix techniques. The results are of in- 
terest both for cosmic strings and for their condensed matter analogues, superfluid 
and superconductor vortices. We tabulate the instabilities of higher winding num- 
ber vortices, and find several bound states. These localised coherent oscillations of 
the order parameter can be thought of as particle states trapped in the core of the 
string. 
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Introduction 



The ideas of gauge unification and spontaneous symmetry breaking lie at the heart of 
modern particle physics. The well established view is of some Grand Unified Theory 
with gauge symmetry group G, suffering a series of spontaneous symmetry breakings 
which reduce the group to that of the Standard Model. In a cosmological context this 
series of spontaneous symmetry breakings is seen to occur as the universe cools. The 
first symmetry breaking occurs at an energy of the order 10 16 Gev, and it is possible that 
topological defects may form, via the Kibble mechanism, at this transition [I], Such 
objects have many important implications: for example, defects of dimension 2 (domain 
walls) and (monopoles) would come to dominate the energy density of the Universe 
soon after a GUT phase transition, which would necessitate a period of inflation if the 
theory admitted such defects. 

When a U(l) symmetry is broken at a phase transition, there arises a network of 
cosmic strings ||, defects of dimension 1, which may provide seeds for large scale structure 
formation in the universe. Because strings are very thin and the cosmological scales of 
interest are many orders of magnitude bigger, they have been approximated by line 
sources with energy-momentum density diag([i, 0, 0, — /i) where fi is the energy density of 
the classical solution. There are, however, a few instances when one is interested in the 
string structure. They include: calculating scattering cross-sections (which are important 
for baryon number violation by strings Q); finding bound states; quantum corrections 
to the energy; and curvature corrections for highly curved string ||. All of these require 
a knowledge of the perturbation spectrum about the string. 

In this paper we find the perturbation spectra in the backgrounds of global and local 
U(l) vortices. Simple U(l) strings provide an ideal test-bed for developing numerical 
techniques for eigenvalues. They are not trivial, like the kink ||, nor are they as geo- 
metrically involved as 3D objects such as the sphaleron |7]]. However, one has to deal 
with all the same problems: gauge choice and numerical stability being the principal 
ones. Furthermore, they are a stepping stone on the way to an efficient technique for 
mapping the stability region of the electroweak string, which we shall describe in a sepa- 
rate publication ||. Our methods are quite straightforward: we find the vortex solution 
using shooting or relaxation algorithms. The numerical profile functions are then used in 
constructing a discretised matrix of second derivatives of the fields. This matrix is then 
diagonalised to find the eigenvalues and eigenf unctions, using the standard linear algebra 
packages incorporated into matlab. 

We find a continuum of scattering states, oscillatory at infinity, for each particle 
species, obeying the expected dispersion relation uj 2 = k 2 + m 2 . There are unstable 
modes for higher winding number strings, corresponding to strings splitting in various 
ways. We tabulate these, and show how different angular momenta in the modes lead 
to different splitting configurations. We also find bound states: N on global strings of 
winding number N (we searched up to iV = 5); but on the gauge strings we found only 
one, which crossed into the continuum for a sufficiently large ratio of the scalar to gauge 
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coupling. The bound states are excitations of the fields at the core of the string, decaying 
exponentially at infinity: one can think of them as particles trapped on the string. 

As is well known, global strings are closely related to superfluid vortices, and local 
ones to flux tubes in superconductors. We conclude by discussing the connection and 
outlining the possible significance of our results for condensed matter physics. 

1 Perturbations about global strings 

The model considered in this section has a complex scalar field with the lagrangian 
density 

£=(^0)*(^0)-A(^-t ? 2 /2) 2 , 

invariant under global U(l) transformations. The field acquires a vacuum expectation 
value, breaking the symmetry and giving rise to a particle spectrum of one Goldstone 
boson and a Higgs boson of mass \/2\rf. 

As is well known, there exist time independent classical vortex solutions to the equa- 
tion of motion 

-<9 2 + 2A(|0| 2 -r/ 2 /2)0 = O. 

The solutions are z-independent and are of the form = (77/ \/2)f(p)e lNd where N is the 
winding number. The energy per unit length of such configurations is given by 

The energy is logarithmically divergent for large p 

E = E core + nrfN 2 In (^j , 

where p c is the core radius. In physical situations there exists a physical cut off, such as 
a container or a nearby anti-vortex, making the energy finite. It can be seen from the 
N 2 ln(R) behaviour that two N = 1 global strings have significantly lower energy than 
one N = 2 global string. The decay of an N = 2 string to two N — 1 strings conserves 
the topological charge and so it would be expected that the decay would occur. Similarly, 
one expects the decay of higher N strings into strings of lower winding number, provided 
the decay conserves the topological charge. Decays which would appear to be possible 
on energetic grounds (eg. one N = 2 string to one N — 1 string) but do not conserve 
the topological charge, do not occur because there is an infinite energy barrier seperating 
the two configurations. So only the N — 1 string is expected to be stable. 

Now consider z- independent perturbations to the classical string background c (x) 
of the form c (a; 1 ) + eS(f)(x l )e~ luJt and 4>*(x l ) + eS(j)*(x l )e luJt where i — 1,2. This gives 
corrections to the energy 

E = E((j) c ) + e 2 E 2 + 0(e 3 ), 
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where 

E 2 = J d 2 x(5(j)*,5(j))M ( 



The matrix M is the perturbation operator, obtained by taking second functional deriva- 
tives of the energy. To find the perturbative modes about the string, which can be 
interpreted as particle states, we solve the coupled eigenvalue problem 



where 

/ -V 2 + 2A(2|0 c | 2 -7? 2 /2) 2A0f 

2A0 2 -V 2 + 2A(2|0 c | 2 -^ 2 /2) 



M = 



The string is a cylindrically symmetric solution, so it is convenient to expand the per- 
turbations in angular momentum states 

<ty = J2 s m e^ N+m ^ 9 and 5(f)* = £ s^e^-^ 6 . 



It is also useful to rescale the coordinates by defining r = ^2\r] 2 p to remove the parame- 
ters A and 77 from the equations. So rescaling and substituting in the string solution and 
perturbation expansion gives the coupled eigenvalue equations for the functions s m and 

^—rr, 



( v -V 2 + ^^ + ^(2 / 2 -l))^ m + i / 2 Sm = cV m 

(-V 2 + ^^ + i(2 / 2 -l)) Sm + i/V m = u 2 s m , 

where u 2 is now dimensionless. To recover the dimensionful values, one should mul- 
tiply by m 2 H = 2Xr] 2 . Note that the eigenvalue equations are the same after complex 
conjugation and changing to — > —to, so we need consider only to > 0. 

If we set (s m ,s*_ m ) to real and imaginary parts, we can seperate the above eigenprob- 
lem into two separate eigenproblems with explicitly real fields 



D 1 + U2f 2 -1) \f 2 



s' \ 9 / s 



2\*J ~ ±) 2J »m _ ,,,2 

1/2 n 1 1/0^2 



j 2 d 2 + |(2/ 2 -i) y v s i m y \s 



(1) 



= ^ T (2) 



V i/ 2 D 2 + |(2/ 2 - 1) ) \ -si m ) " \ -sl m 

where D x = -<9 2 + (to - Nf/r 2 and L> 2 = -d 2 + (to + N) 2 /r 2 . 
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First consider the vacuum where f = 1 and N = 0. Then we have D 1 = D 2 and we 
can diagonalise M with the with the eigenvectors (1, l) T s rn and (1, — l) T s m to give 




(3) 



(4) 



We can see that for the vacuum the solutions are Bessel's functions J m {kr) where for 
(3) k is given by the continuum u 2 = k 2 + 1 and for (4) by u 2 = k 2 . The eigenvalue 
spectrum of (3) corresponds to the spectrum a Higgs particle with rest mass 1 in units 
of 2Xr) 2 , whereas the eigenvalue spectrum from (4) has a continuum down to u 2 = 
corresponding to Goldstone modes. 

Numerical method 

The eigenproblem above was solved with the boundary conditions s m ,s_ m — > as r — > oo 
for various values of m. The profile for a string of winding N was found by solving the 
equation for / using a shooting method. This profile was then substituted into the 
discretized perturbation operator M and the eigenvalue problem was then solved by the 
standard matrix methods used by matlab. This proceedure was performed for linearly 
discretized lattices with 64, 128 and 256 points. The results for the 256 point lattice 
are shown in Tables 1 and 2. The numerical values for the zero modes give an estimate 
of the accuracy of the method. The eigenvalue problem was also solved for a linear 
discretization in the variable p where p = tanh(r). This gives a non-linear discretization 
in r with more points in the core of the string. The results from this are in agreement 
with those quoted in Table 1 to the error given by the values for the zero modes. 

Zero modes 

The string is a solution centered at a point in a two dimensional plane and so there must 
be a perturbed field configuration corresponding to the infinitesimal translation 



For example, for the translation mode in the ^-direction we take the ^-derivative of the 
string solution 



(f>(x) — > (f>(x + Sx) 




<90 c 



cos 



<90 c sin (90c 



dx 



dr r 89 
df iNe iN sin 



= cos6—e 



dr r 



fe 



5 



As can be seen, the translation mode has m — 1, and the forms for s±, s_i can be read 
off from above and shown to satisfy 




The numerical results for this zero mode are shown in the table and, as mentioned above, 
they give an estimate of the error in the eigenvalues. 

Bound states 

For the case m = 0, we have D\ = D2 and again we can diagonalise M with the 
eigenvectors (1, l) T s and (1, — l) T s to give 

(-V r 2 + ^ + i(3/ 2 -l)) So = u 2 s (5) 

(-V 2 + ^ + i(/ 2 -l)) So = (6) 

The eigenvalue spectra of the above equations have continua as for the vacuum, plus 
the appearance of descrete eigenvalues with uj 2 < 1 for (5) due to the potential confining 
states on the string. These later states can be interpreted as Higgs particles trapped on 
the string or as oscillations of the string thickness. The continuum states are scattering 
states with solutions tending to phase shifted Bessel's functions as r — > 00. We found 
that there N bound states for a string of winding number N for iV = 1, ... ,5. The 
solutions to (6) consist solely of scattering states. 

If we consider the perturbations as a function of z then the perturbation expansion is 

5(f) = J2 s m e l{N+m)e e ik * z and 5(f)* = ]T s*_ rn e~ %{N - m)e e~ ik ^ . 

and the solutions to the corresponding eigenvalue problem are similar to those given above 
but with uj 2 = k 2 + k 2 + 1 etc. The bound modes given above can then be interpreted as 
travelling modes along the string. The profile for the bound mode of the N = 1 vortex 
is shown in Figure 1 together with the string profile. 

Decay modes 

As noted earlier all global strings with N > 2 are expected to decay to strings of lower 
winding number. The modes corresponding to these decays occur for m > 2. The 
eigenvalues are listed in Table 1, and we display + e5(f> (with a large e) in Figures 2-5, 
to give a graphical representaion of the decay. 
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For N = 2 there is an m = 2 mode with negative eigenvalue, corresponding to the 
string splitting into two N — 1 strings (Figure 2). 

For N = 3 there is an m = 2 mode with negative eigenvalue corresponding to the 
string splitting into three N — 1 strings in a line (Figure 3). The m = 3 negative mode 
corresponds to the string splitting into three N — 1 strings in a triangular configuration 
(Figure 4). There is also an m = 4 negative mode which gives rise to five vortex cores 
(Figure 5). To see what is happening here consider the small r behaviour of the perturbed 
string solution 

= C + 5(f) 

= fe m + s 4 e m + s_ 4 e~ ie , 

which for small r is given by 

so can see that — > re"* 61 as r — > 0. There is therefore the beginnings of an anti-vortex 
at the origin. Thus the net winding number is three, which is just as well since it is a 
conserved topological charge. 

This is the first example of a more general feature that strings of winding number 
N possess negative modes for 2 < m < 2N — 2 where the string splits into strings of 
lower winding number, m + 1 of them for m > N and for m = N the string splits into 
iV strings of unit winding. If we consider the energy of the string and the conservation 
of topological charge, the general rule for the allowed decays is that an iV wound string 
— > M strings of winding Mj , such that 

E M * = N (7) 

i 

< N 2 (8) 

i 

The first condition is the conservation of topological charge. The second condition comes 
from the iV 2 ln(R) dependence of the energy for large R. The J2i = N 2 possibility is 
ruled out due to the dependence of the core energy on N. If the energy of the string is 
divided up as 



E G oldstone = H^N 2 If 2 

J r 

Em 99S = ^7 rrfr f(l) 2 + I (/2 - 1)2 



then Em ggs is finite whereas Ecoidstone has the ln(i?) dependence for large R. As can be 
seen from Table 1, the sum of Eni gg s for the M strings is greater than that of for the 
string of winding N. 
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Although there are decay modes which begin to generate anti-vortices at the core for 
small perturbations, do they remain as the string decays further? In any physical situa- 
tion this is unlikely. This is because the decay modes with anti-vortices have a number 
of unit winding number vortices arranged arround them in a symmetrical arrangement 
(for instance the m = 6 decay mode of the N = 5 string has a unit winding number 
anti- vortex at the center and six unit winding number vortices arranged arround it in a 
hexagon), and any slight asymmetry will result in the anti- vortex being attracted to one 
of the vortices and annihilating. 

2 Perturbations about local strings 

In this section we are considering the abelian-Higgs model with the lagragian density 
C = - i V" + (D^TiDy) - A(|0| 2 - rf/2f, 

where 

= (d^ + ieA^cj) 
Fftu = d^A u — d u A^. 

The lagrangian is invariant under local U(l) transformations. This symmetry is broken by 
the Higgs field acquiring a vacuum expectation value, and this gives a particle spectrum 
of a massive gauge boson of mass er] and a massive Higgs boson of mass \f2\rf. As 
shown by Nielsen and Olesen || there exist time independent classical vortex solutions 
to the time independent field equations 

D a D a cj> = 2A(|0| 2 - r7 2 /2)0 
d a F ab = -ie((f>*D a ((>-<l>(D a <f>)*), 

where a = 1,2,3. The cylindrically symmetric z-independent vortex solutions can be 
written in the form <fi = (r// y/2)f(p)e lNe , Aq = —Na(p)/ep and A p = where N is the 
winding number of the string. The energy per unit length of the configuration is given 

by 

where the new dimensionless coordinates r is defined by 

erj rriA 

and the parameter (3 is 

m| = 2A 

m A e z 



8 



It will also be convenient to rescale the fields by making the replacements 

Dimensionless fields and coordinates will be used throughout this section. The stability 
of iV > 1 Nielsen- Olesen vortices to splitting into vortices of lower winding has been 
studied previously [[[0], where it was shown that they are stable for j3 < 1 and unstable 



for (3 > 1. This can be seen from Figure 6 showing the energy per vortex for three values 
of f3. For /3 > 1 the energy per vortex is lowest for N — 1, and so a vortex with winding 
number greater than one will split into vortices of lower winding number and so move 
left on the curve. For (3 < 1 the reverse is true. 

Perturbations and gauge fixing 

Now, consider perturbations to the classical string solution of the form (p c (x) +eS(p(x)e~' luJt , 
4>*,(x) + e5(p*(x)e' luJt and A^(x) + e5A fl (x)e~ tuJt . This gives corrections to the action 

S = S(<P c ,A c e )+e 2 S 2 + 0(e 3 ), 

where 

S 2 = ^J d 4 x5&V5$, 

and where T> is the perturbation operator and 5& = (<5</>*e la; ', 5(pe~ tuJt , 5v4 M e~ lw *). The 
perturbation operator T> is given by 

V = 

-pep* 2 -d 2 + 2iA^ + Al - (3{2\(j)\ 2 - I) 2id^* + i(p*d^ + 2A^ 
V id lt <P*-i<P*d ll + 2A fl <ff -id fl <f> + i<pd ll + 2A fl <f> g^d 2 -d^d u + 2g^\(J)\ 2 

and the equations of motion for the perturbations are D5$ = 0. 

In order to set up an eigenvalue problem for the perturbations we need to remove the 
linear derivative terms, as well as removing the gauge degrees of freedom. To do this we 
choose the background gauge condition fllf 



F(A) = d^SA* - (50*0 C - ^M) = 0. (9) 
which is imposed by adding the gauge fixing term 

C GF = \\F{A)\ 2 

to the lagrangian. This enables us to seperate out the time derivatives to give the 
eigenvalue equations 



M 



GF 



1 s * \ 


( W \ 


5(f)* = u 2 




V SAi J 
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(-V + 2W')(£)=<^(£) 

1, 2 and the gauge fixed perturbation operator is 
M GF = 

/ -V 2 - 2^V, + A\ + /3(2|0| 2 - 1) + |0| 2 (0 - 1)0 2 - 2iV l <j) + 2^0 \ 

(/3-l)0* 2 -V 2 + 2^V l + 4 2 + /3(2|0| 2 -l) + |0| 2 2iVi0* + 2A^* ■ 
V 2iVi0* + 2A0* -2iVi0 + 2Ai0 -5 ij V 2 + 25 ij \<f)\ 2 / 

The 5A and 5A Z perturbations decouple because the background string solution is 
independent of t and z. 

The gauge condition (9) above does not fix the gauge completely. For comparison, 
consider the more familiar coulomb gauge V ■ A = ,where it is still possible to make 
gauge transformations 

Ai—>Ai- V;A 

and the gauge transformed field still satisfy the gauge, as long as A is harmonic. Similarly 
for the background gauge condition it is still possible to make gauge transformations for 
A satisfying 

(-«9 2 + 2|0| 2 )A = O. 

This is why in quantum field theory the above gauge fixing term would be accompanied 
by a Fadeev-Popov ghost term of the form 

Cfpg = ^A*(-d 2 + 2|0 c | 2 )A. 

In 2+1 dimensions the perturbation operator for 5A and the ghosts is the same, 
so 5A and A will have the same eigenvalue spectrum. In any calculation of quantum 
corrections, the perturbation operator for the ghosts enters with the opposite sign and a 
factor of 2 compared to that for 5A . Thus the eigenvalue spectrum of the ghosts cancels 
that from 5A and a subset of the spectrum from M GF . We can therefore ignore 8A , 
and we must beware that not all the eigenvalues of M GF correspond to physical states. 

In 3+1 dimensions, where the modes can have a z dependence, the ghosts cancel the 
contributions to the energy from a linear combinations of SAq and SA Z , and the spectrum 
built on the unphysical eigenmodes of M GF . 



Numerical method 

We expand the scalar field perturbations in angular momentum states and the gauge 
field in total angular momentum states given by 

a - A o, , A e td . . ,Aq 

A + = —=(A r -i — ) and A_ = -^(A r + i — ) 
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The background string solution is 



C = fe im , A\ = 



V2r 



A c _ = 



-iNa 

7= ( 

V2r 



and the perturbations are 



i(N+m)e 



50* = s^ m e~ i{N - m)e 



5A+ = ^ia m e 



i(m-l)d 



6A. 



— E * a -; 



Substituting in the above gives the eigenvalue equations 



where 



( D 1 


V 


A 


B 


\ 


( 




\ 




( 




\ 


V 


D 2 


B 


A 






S* 

°-m 




= uj 2 




S* 




A 


B 


Dz 









d m 










V B 


A 







J 


V 


a -m 


) 




\ 




) 



D 1 = 


^2 , (N(l-a)+m) 2 


D 2 = 


^2 , (N(l-a)-m) 2 


D 3 = 


V 2 + (m - 1)2 +2f 


D 4 = 


-V 2 + (m+ 2 1)2 +2f 


V = 


(/3-1)/ 2 


A = 




B = 





+ /?(2/ 2 -l) + / 2 
+ (3(2f 2 -l) + f 2 



As before, the eigenvalue equations are the same after complex conjugation and changing 
m — > — m so we need consider only m > 0. 

If we resolve the functions s m , s*_ m , a m , a*_ m into real and imaginary parts then, as 
for the global string, the complex eigenproblem separates into two eigenproblems with 
explicitly real fuctions 



( D x V A B \ 

V D 2 B A 

A B D 3 

V B A D 4 J 



^ S rn \ 

v 

V 

°— m 

v 
a m 

V a r -m J 



LO 



s r 

" — m 



(10) 
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( Dl 


V 


A 


B 




I 


°m 


\ 






( 




\ 


V 


D 2 


B 


A 










u 2 








A 


B 


D 3 









a m 








a 1 




V B 


A 





D 4 


) 


V 




J 






\ 


~ a -m 


I 



(11) 



where D a , V , A and B are given above. 

First, consider the vacuum, where / = 1, N — and Ai = 0. Then A — B — 0, Di — 
D 2 in the above, and the upper left block can be diagonalised with the eigenvectors 
(1, l) T s m and (1, — l) T s rn to give the four Bessel's equations 

uJ 2 s m (12) 

u 2 s m (13) 

uJ 2 a m (14) 

u 2 a_ m . (15) 

The eigenvalue spectra of these are u 2 = k 2 + 2/3 for (12) corresponding to that of a 
Higgs particle, and oj 2 = k 2 + 2 for the other three corresponding respectively to the 
longitudinal and transverse states of the massive gauge boson. 

In the presence of the string background we still obtain the continua given above, 
but as for the global string, there are discrete eigenvalues correponding to particle states 
trapped on the string. The profiles f(r) and a(r) were solved for with a relaxation 
method on the energy functional, and then substituted into the eigenvalue equations 
above. The eigenproblem was solved with the boundary conditions s m , s_ m , a m , a_ m — > 
as r — > oo for various values of m, by the same method used for the global string. This 
was performed for linearly descretized lattices with 50 and 100 points. The results given 
in Table 3 are for 100 points and the translation modes give an estimate of the accuracy. 
Note that the error in the translation modes is greater for f3 = 4 than f3 — 1 or 0.25. This 
is because the Higgs core is narrower with a much sharper transition from to 1 than 
in the other cases, and so there are fewer points describing this region. One might think 
that a similar situation would hold for the gauge field for (3 = 0.25 but its transition from 
to 1 is more gentle. 

Zero modes 

After gauge fixing the only mode which can have zero eigenvalue is the translation mode. 
The field configurations corresponding to a gauge fixed infinitesimal translation are 

4>(x) -> (f)(x + 5x) = 4>{x) + (D j (f))5x j 



-V 2 + ^ + 2i).s = 



m 



- V" + — + 2 s. 



-V 2 + 



[m - I f 



+ 2) a, 



. (m + 1) 2 \ 
-V 2 + ^ + 2 a_, 
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= (j) + e5(j) 
A, t (x) -> A, t (x + 5x) = Ai(x) + Fji5x j 

= Ai + edAi. 

These perturbations satisfy the gauge condition (9) by the equations of motion. They 
also satisfy 

M5$ = 0. 

by the equations of motion. The numerical results for the translation eigenvalues of M GF 
are given in Table 3. 

Bound states 

Let us first consider m — 0. In this case D\ = D2 and D% = D4, and M has eigenvectors 
of the form (so, So, ^o, a o) T , (so, — $o, — clq) t . Hence there are two sets of solutions for 
so and ao, each satisfying one of two equations: 

D, + V 
A + B 

D 1 -V 
A-B 

For each N, we find that there exists a maximum value of (3 below which equation flT6|) 
possesses one bound state solution, corresponding to some sort of coupled Higgs-gauge 
boson trapped on the string. For the N = 1 string the bound state exists for (3 < 1.5. 
For f3 > 1.5 this bound mode becomes one of a number of modes in the continuum, for 
which the Higgs field is bound but the gauge field is not. These modes correspond to 
some sort of resonant scattering state of the gauge field. For (3 < 1 there are similar 
resonant scattering modes, but with the Higgs and gauge fields interchanged. For higher 
values of N there is still only one bound mode, but the upper limit on j3 is higher. There 
are also several resonant scattering states. 

The solutions to equation (|T7|) also possess bound states for a range of (3. For the 
range 1 < (3 < 8 there is one such mode. For [3 < I, and for higher N, where there 
are a number of bound states. By examining the asymptotic form of the perturbation 
operator, one can establish that the mass squared is 2, and thus we are dealing with 
modes of the gauge field. The gauge mode is trapped on the string because it has a lower 
mass inside the Higgs core, so we would expect the bound mode to have lower energy and 
for there to be more of them when the Higgs core is wider than the gauge field core. This 
is indeed the case. However, as we warned earlier, gauge field modes are not necessarily 
physical. In fact, we find that these vector bound states have very similar eigenvalues 
to the ghost modes, strongly indicating that they are just a result of the residual gauge 
invariance in the background gauge. 




to 



to 




(16) 
(17) 
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Table 3 shows the eigenvalues of the lowest bound modes occuring for equation (|T6|), 
equation ([[7]) and for the ghosts. One can see that the ghost bound modes eigenvalues 
are very close to those from equation (|T7]). The zero mode eigenvalues give an indication 
of the error: the difference between the gauge field and the ghost eigenvalues is actually 
greater than the departure of the numerical zero mode eigenvalue, but we believe that 
our identification of them is nevertheless correct. 

In 3+1 dimensions, there is also the A z field to consider. This has exactly the same 
perturbation equation as the ghosts and SA . As there are in total 2 ghost degrees of 
freedom, they can only cancel one linear combination of 5Aq and SA Z , leaving the other 
for real physical propagating modes, at least one of which is a bound state for (3 < 8. 
Thus we have a massive gauge boson confined to the string. 

Decay modes 

Abelian-Higgs strings with N > 1 for (3 > 1 are unstable to decay to strings of lower 
winding number. The gauge field present in the local strings has the effect of removing 
the energy's logarithmic dependence on the radius, and so the allowed decays of local 
strings are N wound string — > M strings of winding Mj with I]i M; = iV and Mj > 1. 
For (3 < 1 the modes are the same but the eigenvalues become positive. Consider the 
eigenvalues given in Table 3 for the m = 2 mode for N = 2. When (3 = 4 the string 
is unstable to splitting into two N = 1 strings. As is well known (3=1 strings do not 
interact and so uj 2 = for this case. Now when (3 = 0.25 this splitting mode has a small 
positive eigenvalue. This mode then corresponds to two N = 1 strings oscillating about 
their center of mass position. Again if we consider the z direction we obtain travelling 
waves along the string. Table 3 gives the eigenvalues of the splitting modes for various 
(3. 

Discussion 

For global strings of unit winding number, we have found a bound mode with zero 
angular momentum which, when quantised, corresponds to a bound Higgs particle. The 
interaction terms in the lagrangian, however, give a vertex for the Higgs boson to decay to 
two Goldstone bosons, with a decay rate Y ~ m#. A Higgs boson occurs as a resonance 
in the scattering of two Goldstone bosons at a centre of mass energy m 2 H . Similarly, the 
Higgs boson trapped on the string will be seen as a resonance at a centre of mass energy 
of approximately 0.81m 2 H . Indeed, in a strongly coupled theory, the state could only be 
seen as a very broad resonance, rather like the a in tctc scattering. 

For local strings we have found a bound coupled Higgs-gauge boson for (3 < 1.5 and 
a bound gauge boson in 3+1 dimensions for 1 < (3 < 8. There are vertices for a Higgs 
boson to decay to gauge bosons, but for the bound modes this cannot occur due to 
energy conservation. So for U(l) cosmic strings, there are bound modes for the above (3 
ranges. Populations of these bound states could have interesting effects: if they are truly 
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stable, they could alter the equation of state of the string, adding to the mass density 
and subtracting from the tension. However, in a realistic model it is rather unlikely that 
the string fields will not be coupled to light fields, which offer a decay route to the bound 
states. 

Global strings possess a condensed matter analogue in the Ginzburg-Pitaevski theory 



12, Oj of superfluid Helium II. In the G-P theory the free energy of the condensate is 



given by 



d 3 x &VVf -«h/f + fl 



where if) is the condensate wave function, which for stationary superfluid has \ipo\ 2 = 
Ps/ m i — ac/0, where p s is the superfluid density and m 4 is the mass of a helium atom. 
The equation of motion is first order in time, commensurate with the non-relativistic 
nature of the system: 

dip _ SF 

There exists an extremum of the free energy of the form f(r)e lNe where f(r) is the same 
as in Section 1. This solution decribes a vortex in the superfluid in which the superfluid 
flows at velocity v s = UN/ m^r about a line along the z-axis, where the superfluid density 
is zero. This is analogous to the global cosmic string. Unfortunately, we cannot rely too 
heavily on this description of vortices in Helium II, because the core size is of the same 
order as the inter-atomic spacing - the theory is strongly coupled. It is nevertheless 
interesting to expand about the vortex solution to find the spectrum of states in the 
vortex background. To make contact with the dimensionless version of the perturbation 
matrix M of Section 2, it is convenient to scale all distances in units of the correlation 
length £ = v /(^ 2 /2m4«), and to use a dimensionless field x — ^ViP/®)- Then for 
perturbations 5\ we find 





= «(|^o|T)M 

Thus the energy eigenvalues in the matrix background of states with momentum hk z /^ 
along the string are 

e = a(\ip \ 2 e)^ 2 + kl). 

Of particular interest is the spectrum of bound states at u 2 ~ 0.81, which we believe 
has not been convinvingly established before in Helium II vortices. These are not the 



"bound" states found by Pitaevski [JL4]], which are modes with m = 1 associated with 
the translation of the vortex as a whole. What we are discussing is a zero angular 
momentum oscillation of the size of the vortex - a "sausage" mode. Bound states at 
angular momentum and 2 have been reported by Grant [fTS!, but the accuracy of the 
numerics is open to question, for both states were found very close to threshold (u 2 = 1). 
We find no evidence of bound states at m — 2, and at m = u 2 is approximately 20% 
below threshold. 
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These bound modes may be observable as a resonance in the scattering of phonons 
or second sound off the core of the vortex, but in the light of the comments above about 
the applicability of G-P theory to such a vortex, this may not be observable. 

The Abelian-Higgs model also has a condensed matter analogue in the Ginzburg- 
Landau theory |IB| of superconductivity. It has been shown that there is a range of 



conditions for which G-L theory is applicable JT7 
In G-L theory the free energy is given by 

where ip is condensate wave function and B is the magnetic field. For the condensate 
in the absence of magnetic fields \ipo\ 2 = a/ (3. There are two length scales to use in 
constructing a dimensionless free energy functional, the correlation length £, and the 
penetration depth A = v /(/3m/4e 2 /i «)- The Ginzburg-Landau parameter is defined as 
the ratio k = A/£. Defining a dimensionless gauge field V{ = 2eAAi/h, we find 

F = ^«(|^ | 2 A 3 ) J d 3 y (|(V„ - iv a ) X \ 2 + ^(V x vf - 2k 2 \ X \ 2 + k 2 \ X \ 

We see that the Ginzburg-Landau parameter is analogous to i/(/?/2), where (5 is the ratio 
of scalar to gauge couplings used earlier. This gives the same perturbation matrix M 
(and M ), up to an overall scale factor, and so the energy eigenvalues of the states in 
the superconductor vortex background are 

- -^-M\M 2 ^ 2 + k 2 z ). 



2k 

The bound modes on the string exist for a range of k. Converting the j3 limits into n 
gives that for l/\/2 < k < a/3/2 there are bound modes of type (a) (coupled Higgs- vector 
modes), while for \j\[2 < k < 2 there are bound modes of type (b) (longitudinal gauge 
field modes). 

Again, the observability of such modes is an open question, complicated by the fact 
that superconductor vortices are already known to have quasi-particle bound states [IT 
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Table captions 



Table 1: Emggs of global strings with winding number N in units ixrf and m = bound 
mode eigenvalues in units of 2\rf. 

Table 2: Eigenvalues of decay modes of global strings of winding number N, in units 
of 2\rf. m is the angular momentum of the mode and the Mj values are the winding 
numbers of the strings formed by the decay. 

Table 3: Eigenvalues for local string modes in units of e 2 i] 2 /2. The modes shown are 
a)Higgs-gauge boson bound mode, b)gauge boson bound mode, c)ghost bound mode 
t)translation mode, and s)splitting modes. 

Figure captions 

Figure 1: Global string field profile f(r) for iV = 1 and the m = bound mode field 
profile Sq(t). 

Figure 2: m = 2 decay mode for a N = 2 global string. 
Figure 3: m = 2 decay mode for a N = 3 global string. 
Figure 4: m = 3 decay mode for a N = 3 global string. 
Figure 5: m = 4 decay mode for a N = 3 global string. 

Figure 6: The energy per vortex for N — 1, . . . , 5 local strings in units of 7irj 2 . 
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N 


^Higgs 


zero mode 


u) 2 of m: 


=0 bound mode 




1 


0.7833 


0.0012 


0.8134 






2 


2.4149 


0.0019 


0.6299 


0.9554 




3 


4.9835 


0.0013 


0.5180 


0.8747 0.9727 




4 


8.5712 


0.0013 


0.4437 


0.7926 0.9314 


0.9822 


5 


13.0310 


0.0014 


0.3907 


0.7206 0.8835 


0.9527 0.9976 



Table 1: 



N 


m 


Mi 




a; 2 


2 


2 


1,1 


2 


-0.2093 




3 


2 


1,1,1 


3 


-0.1283 






3 


1,1,1 


3 


-0.3289 






4 


1,1,1,1,-1 


5 


-0.1207 




4 


2 


1,2,1 


6 


-0.0889 






3 


1,1,1,1 


4 


-0.2429 






4 


1,1,1,1 


4 


-0.3889 


-0.0617 




5 


1,1,1,1,1,-1 


6 


-0.2410 






6 


1,1,1,1,1,1,-2 


10 


-0.0749 




5 


2 


1,3,1 


11 


-0.0662 






3 


1,1,1,2 


7 


-0.1908 






4 


1,1,1,1,1 


5 


-0.3151 


-0.0369 




5 


1,1,1,1,1 


5 


-0.4224 


-0.1643 




6 


1,1,1,1,1,1,-1 


7 


-0.3145 


-0.0272 




7 


1,1,1,1,1,1,1,-2 


11 


-0.1868 






8 


1,1,1,1,1,1,1,1,-3 


17 


-0.0475 





Table 2: 
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N 


m 


mode 


UJ 2 








P = 4 


P = l 


(3 = 0.25 


1 





a 




1.5540 


0.4538 






b 


1.9097 


1.5302 


1.1073 






c 


1.9304 


1.5565 


1.1298 




1 


t 


-0.0459 


-0.0187 


-0.0096 


2 





a 




1.0766 


0.3604 






b 


1.5110 


1.0476 


0.6775 






c 


1.5713 


1.0781 


0.6963 




1 


t 


-0.0267 


-0.0050 


-0.0001 




2 


s 


-0.7654 


-0.0007 


0.1121 


3 





a 


1.7085 


0.8043 


0.2941 






b 


1.2163 


0.7856 


0.4860 






c 


1.2720 


0.8057 


0.4932 




1 


t 


-0.0152 


-0.0013 


0.0016 




2 


s 


-0.4059 


-0.0009 


0.0617 




3 


s 


-1.3497 


0.0023 


0.1952 



Table 3: 
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Figure 2. 




Figure 3. 




Figure 4. 




Figure 5. 




